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ABSTRACT 


The  charge  distribution  of  periodic  electron  beam  pulses  generates  patterns 
of  sub-Cerenkov  radiation  distinctive  of  the  distribution  of  charge  with  in  a 
bunch.  Mapping  the  radiation  pattern  from  different  charge  shapes  may  provide 
insight  into  whether  charge  pulse  shapes  can  be  determined  from  observed 
radiation  patterns.  The  radiation  patterns  of  Gaussian,  Level,  and  Trapezoidal 
function  were  mapped  by  computer  simulation  Near  90°  to  the  beam,  the 
radiation  patterns  of  all  three  charge  distributions  developed  an  envelope 
proportional  to  the  fourier  transform  of  the  charge  bunch  distribution  when  the 
wavelength  of  the  emitted  radiation  was  comparable  to  the  size  of  the  bunch. 
For  the  Gaussian  function,  the  envelope  is  Gaussian,  for  the  level  function  it 
is  a  sine  function.  Since  the  envelope  for  the  trapezoidal  function  in  the 
product  of  two  sine  functions  it  is  more  difficult  to  analyze.  This  work  may 
provide  a  basis  for  determining  the  charge  shape  of  electron  beam  pulses  from 
the  sub-Cerenkov  radiation  based  on  the  radiation  intensity  pattern. 
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I.  INTRODUCTION 


A.  HISTORY 

A  very  faint  emission  of  bluish-white  light,  what  was  later-on  called 
Cerenkov  radiation,  from  transparent  substances  placed  nearby  a  radioactive 
sources  had  been  observed  by  Mme  Curie  in  1910  .  Although  electromagnetic 
theory  had  been  sufficiently  developed  at  that  time  to  describe  this 
phenomenon,  the  first  deliberate  attempt  to  study  the  phenomenon  was  made 
by  Mallet  after  many  years,  who  found  that  the  light  emitted  from  a  wide 
variety  of  transparent  bodies  placed  close  to  radioactive  source  always  had  the 
same  bluish-white  light  quality,  and  that  spectrum  was  continuous,  not 
possessing  the  line  or  band  structure  characteristic  of  fluorescence  [Ref.  1:  p. 
!]• 

P.A.  Cerenkov  carried  out  a  series  of  experiments  between  1934  and  1938 
related  to  the  phenomenon.  In  1937  Ilya  Frank  and  Igor  Tamm  proposed  a 
satisfactory  theory  of  the  radiation.  The  experimental  results  and  theoretical 
predictions  were  in  excellent  agreement.  A  complete  description  of  C°renkov 
radiation,  as  well  as  exact  mathematical  treatment  is  given  in  [Ref.  1:  p.  1[. 

B.  WHAT  IS  THE  CERENKOV  RADIATION  ? 

A  charged  particle  or  a  group  of  charged  particles  moving  through  a 
medium  at  a  speed  greater  than  the  (phase)  velocity  of  light  in  the  same 
medium  will  generate  a  continuous  spectrum  of  electromagnetic  radiation.  This 
process  is  referred  to  as  Cerenkov  radiation. [Ref.  2:  p.9) 


C.  PREVIOUS  RESEARCH  AT  N.P.S. 

The  research  on  Cerenkov  radiation  has  been  conducted  by  Professors  Fred 
R.  Buskirk,  John  R.  Neighbours  and  others  in  the  Naval  Postgraduate  school 
since  the  LINAC  was  established.  Recent  work  has  concentrated  on  the 
determination  of  electron  beam  bunch  profile  through  Cerenkov  radiation. 

The  work  was  aimed  to  illustrate  examples  of  the  consequences  of  the 
finite  spatial  charge  distribution  to  show  that  Cerenkov  radiation  emitted  at 
low  microwave  frequencies  may  be  used  to  characterize  the  properties  of  an 
intense  relativistic  beam  in  air.  The  results  were  that  the  spatial  charge 
distribution  of  an  electron  pulse,  along  with  the  beam  interaction  length, 
determines  the  Cerenkov  radiation  distribution  as  a  function  of  frequency.  An 
angular  distribution  of  the  Cerenkov  radiation  can,  in  principle,  measure  its 
spatial  charge  distribution. [Ref.  3:p.  1994]. 

D.  PURPOSE 

The  charge  distribution  of  a  bunch  influences  the  pattern  of  radiation 
produced.  Detailed  calculations  for  several  simple  charge  distributions  have 
been  carried  out  for  the  Cerenkov  case  [Ref.  4:pp.  1994,1996].  The  object  of 
this  work  is  to  assess  the  effect  of  the  charge  distribution  of  a  bunch  on  the 
sub-Cerenkov  radiation  from  a  charged  particle  beam.  Specifically,  this  means 
investigation  of  the  short  beam  patF  region  at  energies  below  the  threshold. 


II.  THEORY 


A.  BACKGROUND 

A  fast,  charged  particle  may  cause  electromagnetic  radiation  in  a  dielectric 
medium  when  it  moves  with  a  constant  speed.  Suppose  an  electron  to  be 
moving  relatively  slowly  through  a  transparent  medium  as  shown  on  Fig.  1.1 
a).  As  the  particle  passes  through  the  medium,  each  elemental  region  of  the 
medium  along  the  track  will  in  turn  receive  a  very  brief  electromagnetic  pulse 
and  be  polarized,  producing  a  dipole.  An  electromagnetic  pulse  is  generated  by 
the  formation  and  disappearance  of  the  dipole  as  the  bunch  passes.  The 
electromagnetic  pulse  propagates  away  from  the  dipole  source  at  the  velocity  of 
light  in  the  medium.  The  distance  between  wave  fronts  is  compressed  as  the 
velocity  of  the  electron  bunch  increases.  This  effect  is  similar  to  the  doppler 
effect  for  sound  radiated  from  a  moving  source.  If  the  speed  of  the  particle  is 
less  than  the  propagation  velocity  in  the  medium,  the  electromagnetic  field 
produced  is  called  sub-Cerenkov  radiation.  [Ref.  1:  pp.  3-4] 

If  the  electron  is  moving  faster  than  light  in  the  medium  (see  Fig  1.1  b). 
the  the  pulse  source  velocity  is  greater  than  that  of  the  propagating 
electromagnetic  wave  which  transport  the  energy.  The  result  is  a  wavefront 
and  where  the  emanated  wavelets  bunch  together  in  phase.  This  shockfront  of 
electromagnetic  radiation  is  Cerenkov  radiation.  [Ref.  1:  pp.  4-5] 


B.  CHARACTERISTICS  OF  CERENKOV  RADIATION 

Fig.  2.2  shows  the  resultant  Cerenkov  radiation  plane  wavefront  BC  from 
particles  traveling  from  A  to  B.  Wavelets  from  point  0  to  5  are  coherent 

along  wavefront  BC.  The  angle  Qc  is  given  by  the  "Cerenkov  relation": 

cos  =  1  /  3  n,  (2.1) 

where  n  is  the  refractive  index  of  the  medium  and  3  is  the  ratio  of  the 

velocity  of  the  particle  to  the  velocity  of  light  in  a  vacuum  [Ref.  1:  pp.4-5]. 
From  the  above,  the  following  is  seen  that: 

(1)  For  medium  of  given  refractive  index  n.  there  is  a  threshold  velocity 

^min  ~  ( 1/n).  below  which  no  radiation  tak«'  place  [Ref.  1:  p.  5]. 

(2)  For  ultra-relativistic  particles,  for  which  3  i  1. there  is  a  maximum 

angle  of  emission,  given  by  0mAX  =  cos‘'(l/n). 

(3)  Radiation  occurs  mainly  in  the  visible  and  near  visible  region  for 
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which  n  >  1.  In  the  X-ray  region,  above  atomic  resonances,  n  is 
usually  less  than  1;  therefore  Cerenkov  radiation  in  the  X-ray  region 
is  not  probable  [Ref.  3:  p.  3252,  Ref.  1:  p.  6]. 


(4)  In  a  gas.  0C  has  only  a  slight  dependence  on  electron  velocity  because 
,i  must  be  close  to  one  in  order  to  obtain  the  Cerenkov  condition 
[Ref.  3:  p.  3249], 


From  the  nature  of  the  radiation,  which  has  been  described  above,  it  may  be 


concluded: 

(5)  Azimuthal  symmetry  of  the  radiation  produces  a  cone  of  radiation. 

(6)  The  distribution  in  0  of  the  radiation  intensity  approximate  a  Dirac 
^-function  centered  about  Cerenkov  angle  '6C'- 

(7)  Length  AB  in  Fig.2  must  be  much  larger  than  the  radiated 
wavelength  in  order  to  avoid  diffraction  effects.  For  a  finite  radiation 
length,  the  diffraction  effect  is  significant. 

C.  CALCULATED  POWER  OF  COHERENT  CERENKOV  RADIATION 

The  total  coherent  power  per  unit  solid  angle,  radiated  at  the  frequency  v 
by  a  periodic  charged  particle  beam  in  traveling  a  finite  distance  L  at 
constant  velocity  is 

W(i/,k)  =  Q  R2,  (2.2) 

where  vq  is  the  fundamental  frequency  of  the  beam  generator  and  v  is  a 
harmonic  of  i/o,  q  is  the  charge  of  an  individual  bunch,  Q  is  a  constant  and  R 
is  the  radiation  function.  The  bunches  in  the  beam  are  assumed  to  be 
unchanging  in  shape  and  size  as  the  beam  travels  through  the  medium  (usually 
air)  at  a  velocity  given  by  v  =  0  Co-  The  velocity  of  light  in  vacuum  is  Co 
and  the  velocity  of  light  in  the  medium  is  c  =  Co/n  where  n  is  the  index  of 
refraction.  The  parameters  used  to  describe  the  radiation  are 

Q  =  /icq2/8?r2  (2.3) 


14 


(IKi 


♦ 

} 

t 

which  is  a  constant,  where  p  is  the  permeability  of  the  medium,  and  the  J 

i 

radiation  function  R  is  given  by  j 

P 

i 

i 

R  =  2xt]  sin#  I(u)F(k),  (2.4)  ! 

t 

i 

t 

I 

where  6  is  the  angle  between  the  direction  of  travel  of  the  charged  particle  ; 

beam  and  the  direction  of  propagation  of  the  emitted  radiation,  I(u)  is  the  ! 

f 

diffraction  function,  and  F(k)  is  the  form  factor.  [Ref.  5:  pp.2741,2742]  | 

I 

The  wave  vector  of  the  emitted  radiation  is  k  (k  =  u/c),  and  F(k)  is  the 
dimensionless  form  factor.  That  is,  if  p{ r)  is  the  charge  distribution  of  a  single 
bunch,  then  the  Fourier  components  of  the  charge  are 


00 

p(V)  =  |  |  |  p(r)eik'rd3r, 
—  00 


(2.5) 


and  the  form  factor  is  defined  by 

p(  k)  =  qF(k).  (2.6) 

For  a  point  charge,  F(k)  is  identically  one  for  all  values  of  k. 

The  dimensionless  beam  length  parameter  rj  which  is  defined  as  the  ratio  of 
the  path  length  of  the  charged  particle  beam  to  the  wavelength  in  the  medium 
of  the  emitted  radiation,  appears  in  the  radiation  function  explicitly,  and 
implicitly  as  part  of  the  diffraction  parameter  defined  below: 

r]  ~  L  /  A.  (2.7) 
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The  diffraction  function  is 


I(u)  =  sin  u/u,  (2.8) 

where  the  dimensionless  diffraction  parameter  u  depends  upon  the  Cerenkov 
angle  given  by  cos  9C  =  (n/?)"1,  as  well  as  upon  the  beam  length  parameter 
and  the  radiation  emission  angle: 

u  --  7r  7 7  [  (n/3)'1  -  cos  0  ].  (2.9) 

The  parameter  rj  determines  the  periodicity  of  the  diffraction  pattern  and  also 
determines  its  magnitude.  Since  R  depends  directly  upon  rj,  the  strength  of  the 
radiation  is  proportional  to  the  square  of  the  path  length  of  the  beam.  [Ref.  5: 
pp.  2741,2742] 

D.  FINITE-LENGTH  PATH 

If  the  emission  region  has  finite  length,  which  may  be  realized  by  passing 
the  electron  beam  through  a  gas  cell,  the  radiation  propagation  direction  is  not 
confined  to  a  sharp  Cerenkov  angle  9C,  but  is  spread  over  a  range  of  emission 
angles.  For  the  case  of  gas,  for  which  the  ordinary  Cerenkov  angle  is  a  few 
degrees,  the  broadening  is  asymmetric  about  the  Cerenkov  cone,  and  the 
radiated  power  displays  various  interference  lobes.  This  spreading  effect  depends 
only  on  the  path  length  in  the  gas  and  does  not  depend  on  the  electron  bunch 
structure  or  periodicity  of  the  bunches.  [Ref.  3:  pp.  3246,3251] 

Radiation  from  the  polarized  molecules  along  the  path  of  the  electron 
beam  is  formally  equivalent  to  diffraction  from  a  single  slit  at  angles  far  from 
normal  incidence.  The  shift  in  the  radiation  peak  can  be  understood  as  a 


diffraction  effect  arising  from  the  linearly  varying  phase  of  the  radiation  along 
the  beam  interaction  length  [Ref.  3:  p.  3246]. 

Considerable  power  (about  half  in  the  case  of  a  gas)  appears  outside  the 
main  lobe,  and  the  latter  is  the  peaked  at  angles  much  larger  than  the 
Cerenkov  angle.  The  Cerenkov  power  radiated  per  unit  path  length  is  increased 
by  large  factors  approaching  2  orders  of  magnitude  for  a  gas  path  of  finite 
length  compared  to  one  of  infinite  length  [Ref.  3:  p.  3251]. 

This  increase  in  power  is  associated  with  the  finite  length  of  the  radiating 
medium  only.  The  effect  should  be  approximately  the  same  for  either  single  or 
periodic  electron  bunches  and  should  occur  whether  the  bunches  are  effectively 
point  charges  or  have  significant  size. 

The  value  of  the  parameter  u  is  affected  by  77  and  the  beam  energy 
through  0.  As  either  of  these  parameters  is  varied,  the  angular  position  of  the 
zeros  of  I(u)  changes.  The  emission  threshold  is  defined  as  the  situation  when 
the  first  zero  of  the  diffraction  function  enters  the  physical  range.  Applying 
this  criterion,  the  result  is  that  for  small  values  of  77,  a  Cerenkov-like 
radiation  pattern  is  produced  even  when  the  product  n 0  <  1.  That  is,  the 
threshold  occurs  for  n0  <  1  at  small  values  of  77.  As  77  increases  (at  constant 
n0  <  1),  the  threshold  is  crossed  and  the  radiation  pattern  becomes  more 
chaotic.  A  full  discussion  of  these  geometrical  effects  is  given  in  Ref.  5. 

E.  PERIODIC  ELECTRON  BUNCHES 

Cerenkov  radiation  was  considered  for  periodic  bunches  of  electrons  such 
as  would  be  emitted  by  a  typical  traveling-wave  electron  linear  accelerator 
(LINAC).  Because  the  distribution  of  intensity  of  Cerenkov  is  proportional  to 
the  frequency,  the  radiation  at  microwave  frequencies  would  be  low  unless 


j 


beams  are  intense  and  bunched  so  that  coherent  radiation  by  many  elections 
contributes  [Ref.  2:  p.  1531]. 

Periodic  bunches  of  electrons  produce  Cerenkov  radiation  at  harmonics  of 
the  bunch  frequency,  while  a  single  bunch  will  radiate  with  a  continuous 
frequency  distribution.  The  effect  of  bunch  size  is  reflected  in  the  factor  of  the 
Fourier  transform  of  the  charge  in  the  bunch,  which  causes  the  radiation  to 
fall  off  at  high  frequencies,  while  at  low  frequencies  all  charges  in  the  bunch 
radiate  in  phase. 

At  low  frequencies  such  that  the  wavelet  jth  of  the  emitted  radiation  is  of 
the  order  of  the  bunch  size,  the  electrons  in  the  bunch  radiate  coherently.  This 
leads  to  large  enhancement  factors  for  typical  bunches  consisting  of  108 
electrons,  and  in  fact  allows  the  radiation  to  be  significant  at  microwave 
frequencies.  [Ref.  3:  p.  3246] 

At  higher  frequencies  or  harmonics,  destructive  interference,  described  by 
the  Fourier  transform  of  the  charge  density  decreases  intensities  with  increasing 
frequency,  until  incoherent  radiation  takes  over  when  the  wavelength  of  the 

radiation  is  much  less  than  electron  spacing  [Ref.  3:  p.  3246]. 

When  a  single  electron  passes  through  a  dielectric  medium  where  a 
spatially  periodic  term  is  added  to  the  dielectric  constant,  the  result  is 
radiation  occurring  even  for  electrons  which  do  not  exceed  the  velocity  of 
light  in  the  medium  and  at  angles  other  than  the  Cerenkov  cone  angle.  The 
non-Cerenkov  part  of  the  radiation  is  attributed  to  transition  radiation.  If  the 
electron  velocity  were  lower  so  that  v/c  were  close  to  but  less  than  unity,  the 

peak  in  I  would  be  pushed  to  the  left  toward  0°  of  9  such  that  cos  9C  =  v/c 

would  be  larger  than  1.  But  the  tails  of  the  diffraction  function  I  could  extend 
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into  the  physical  range  1  =  cos  9  =  -1.  and  this  would  be  called  transition 
radiation  and  be  ascribed  to  the  passage  of  the  electrons  through  the 
boundaries  of  the  gas  cell.  [Ref.  2:  p.  1535] 

F.  PROFILE  DETERMINATION  OF  A  ELECTRON  BEAM  BUNCH 

1.  Effect  of  The  Form  Factor 

The  bunch  form  factors  are  Fourier  Transforms  F(k)  which  differ 
depending  on  the  charge  distribution  within  a  bunch.  The  form  factor  is 
identically  one  for  a  point  charge,  and  for  a  finite  charge  distribution  F(k)  = 
1  for  k  =  0.  If  all  the  charge  within  a  bunch  is  of  the  same  sign  F(k)  must 
fall  off  as  a  function  of  k  near  the  origin.  Here  we  consider  only  line  charge 
distribution  and  consequently  the  form  factors  are  functions  of  the  product  of 
kz  =  kcos  9  and  appropriate  lengths  which  are  characteristic  of  the  charge 
distribution  of  a  bunch. 

Thus  the  behavior  of  the  form  factor  F  over  the  range  of  emission  angles 
is  dependant  on  the  ratio  of  the  charge  bunch  lengths  to  the  wavelength  of 
the  emitted  radiation.  For  a  small  value  of  these  ratios,  the  value  of  F  is 
essentially  unity  for  all  angles.  Then  the  emitted  radiation  pattern  is 
modulated  only  by  the  geometrical  envelope  function  as  described  in  Ref.  5. 

However  for  larger  values  of  the  ratios,  the  value  of  F  over  the  range  of 
emission  angles  can  be  significantly  less  than  one.  In  this  case  the  dependence 
of  the  form  factor  on  angle  is  dramatic  since  the  cosine  term  goes  to  zero  at 
ninety  degrees  causing  a  form  factor  of  unity.  If,  at  same  time  the  emission 
length  ratio  is  small,  a  condition  could  be  achieved  where  the  majority  of  the 
Cerenkov  power  is  radiated  in  a  direction  normal  to  the  electron  beam  path. 
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A  complete  angular  map  of  Cerenkov  radiation  patterns  may  be  used  to 
characterize  the  properties  of  relativistic  beam  pulses,  and  may  enable  one  to 
determine  the  beam  pulse  charge  distribution  from  measured  radiated  patterns 
[Ref.  4:  p.  1996]. 

2  Bunch  Shapes  and  Corresponding  f  orm  Factors 

Figure  2.3  diagrams  the  beam  pulse  shapes  and  presents  the 
corresponding  form  factor  expressions.  [Ref.  S:  pp.  16,17] 
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Figure  2.3  Electron  Beam  Charge  Bunch  Shapes  and 
Corresponding  Form  Factors. 
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Figure  2.3  (cont.)  Electron  Beam  Charge  Bunch  Shapes  and 

Corresponding  Form  Factors. 
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Figure  2.3  (cont.)  Electron  Beam  Charge  Bunch  Shapes  and 

Corresponding  Form  Factors. 
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III.  CALCULATION  AND  RESULTS 


A.  PARAMETERS  FOR  CALCULATION 

Professor  J.  R.  Neighbours  and  others  of  the  Department  of  Physics. 
Naval  Postgraduate  School,  Monterey,  California,  developed  a  computer 
program  to  calculate  and  plot  angular  maps  of  the  radiated  energy  per  unit 
solid  angle  for  a  given  form  factor.  The  program  originally  was  designed  to 
run  on  a  Tektronix  Computer  and  considered  the  following  one  dimensional 
(i.e.  in  Z  only)  beam  charge  distributions:  Gaussian,  level  (i.e.  boxcar), level 
plus  a  sinusoidal,  ripple  combination, and  a  double  hump. 

This  work  was  aimed  to  modify  the  program  to  run  on  the  IBM  370 
mainframe  computer  and  plot  the  calculated  power  per  solid  angle,  radiated  at 
the  frequency  i/,  using  the  DISSPLA  V10.5.  The  results  are  limited  to  the 
three  form  factors  for  the  rectangular,  Gaussian  and  trapezoidal  charge 
distributions.  Calculations  for  other  charge  distributions  are  left  for  future 
work. 

To  calculate  the  total  coherent  power  per  unit  solid  angle,  radiated  at  the 
frequency  v  by  a  periodic  charged  particle  beam  in  traveling  a  finite  distance 
L  at  constant  velocity,  a  set  of  parameters  were  chosen  close  to  those 
available  in  an  experiment. 

(1)  Index  of  refraction,  n  =  1.000268. 

(2)  Speed  of  light  in  vacuum.  Co  =  2.997925  *  108  m. 

(3)  Fundamental  frequency,  i'o  =  2.85  GHz. 

(4)  Bunch  charge,  B  =  1.5  *  10'12  Coulomb. 


r.1  u"^y>.'  ^iv-\r’\r,  v. 


These  values  give  an  nd  product  of  n 3  =  0.98,  which  is  2  %  below  the 
infinite  path  length  threshold  value  of  unity.  By  choosing  different  radiation 
harmonics  and  varying  the  path  length  the  radiation  pattern  above  (  i]  =  12  ) 
and  below  (  77  =  37  )  the  threshold  for  2.555  Mev  electron  bunches,  could  be 
investigated.  The  calculations  were  carried  out  at  constant  77  for  the  various 
charge  distribution  as  the  size  of  the  distribution  was  varied. 


B.  GAUSSIAN  FUNCTION 

Fig.  3.1  through  Fig.  3.8  show  the  radiation  intensity  pattern  of  Gaussian 
line  charge  distributions  in  sub-Cerenkov  radiation  (nd  =  0.98  )  for  two 
different  value  of  77.  Fig.  3.1  through  Fig.  3.4  are  the  results  of  using  77  =  12. 
with  N  =  3,  varying  the  bunch  length  from  0.24  cm  to  2.4  cm.  The  rest  of 
them  are  for  the  case  of  using  77  =  37.  and  likewise,  varying  the  bunch  length. 
BL  is  the  symbol  used  for  bunchlength. 

Figures  3. 1-3.4  show  the  effect  on  the  emitted  radiation  pattern  as  the 
bunch  length  is  increased  while  the  beam  path  parameter  77  is  held  constant  at 
77  =  12.  With  a  bunch  length  of  0.24  cm  the  radiation  pattern  is  essentially 
that  of  a  periodic  set  of  point  charges.  As  the  bunch  length  increases,  a  hump 
first  becomes  evident  in  Fig.  3.2  when  the  wavelength  of  the  emitted  radiation 
begins  to  become  comparable  to  the  size  of  the  bunch.  Further  increase  in 
bunch  size  shows  an  increasing  dominance  of  the  form  factor  until  finally  the 
radiation  pattern  is  almost  completely  dependant  on  the  form  factor.  In  Fig. 
3.4  the  envelope  of  the  radiation  pattern  is  the  form  factor  for  the  Gaussian 
charge  distribution  except  for  a  slight  skewing  caused  by  the  geometrical 
envelope  function. 
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Fig.  3.3  Gaussian  Function  with  BL/A  =  0.503,  77  =  12 


Fig.  3.4  Gaussian  Function  with  BL/A  =  0.686, 
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Fig.  3.7  Gaussian  Function  wi 
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Fig.  3.8  Gaussian  Function  with  BL/A 
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Similar  remarks  can  be  made  for  the  longer  finite  path  results  shown  in 
Fig.  3.5  -  3.8  for  rj  =  37.  In  this  set  of  figures,  because  this  value  of  rj  causes 
the  radiation  to  be  below  threshold,  the  radiation  patterns  are  wildly 
oscillating  -  almost  chaotic.  However  exactly  the  same  thing  happens  to  the 
envelope  of  the  radiation.  Fig.  3.8,  calculated  for  a  bunch  length  of  2.4  cm.  as 
is  Fig.  3.4,  also  shows  the  radiation  envelope  to  be  the  form  factor  for  the 
Gaussian  charge  distribution. 


C.  LEVEL  FUNCTION 

Fig.  3.9  through  Fig.  3.22  show  the  radiation  intensity  lobe  development 
as  the  level  function  charge  distribution  length  varies  for  the  same  two  values 
of  t].  Fig.  3.9  through  Fig.  3.15  are  calculated  for  rj  =  12.  the  remainder 
through  Fig.  3.22  for  rj  =  37.  LFL  is  the  symbol  for  the  level  function  length. 

In  Fig.  3.9,  where  the  level  function  length  is  much  smaller  than  the 
radiation  length,  the  radiation  pattern  is  almost  the  same  shape  as  Fig.  3.1 
which  is  for  the  Gaussian  charge  distribution.  Proceeding  through  to  Fig.  3.13, 
as  the  length  of  the  charge  distribution  is  increased  the  form  factor  becomes 
prominent.  Fig.  3.13  is  of  interest  in  that  it  also  shows  considerable  radiation 
at  small  angles.  This  result  occurs  because  the  form  factor  for  the  level 
function  has  secondary  maximum  which  at  low  angles  combine  with  the 
relatively  large  value  of  the  diffraction  function  to  give  substantial  radiation. 

Further  increases  in  level  function  length  give  an  increased  role  to  the 
form  factor  until  finally  in  Fig.  3.15,  it  has  become  dominant.  Here  the 
envelope  of  the  radiation  is  seen  to  be  the  (skewed)  form  factor  of  the  level 
function.  The  small  lone  peak  at  approximately  ten  degrees  occurs  at  the 
maximum  of  the  geometrical  envelope  function. 
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Fig.  3.20  Level  Function  with  LFL/A 
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Fig.  3.16  to  Fig.  3.22  drawn  for  /?  =  37  show  the  increasing  dominance  of 
the  form  factor  below  threshold.  The  pattern  oscillates  very  rapidly,  but 

otherwise  the  results  are  as  above.  In  fact  the  rapid  oscillations  make  the  form 
factor  more  evident  in  this  series  of  Figures. 

D.  TRAPEZOIDAL  FUNCTION 

Both  the  Gaussian  and  level  function  distributions  have  a  single 
characteristic  length;  the  trapezoid  is  more  complicated  in  having  two.  the 

length  of  the  base,  lb,  and  the  length  of  the  top,  lt.  In  the  fourier  transform 
these  appear  as  t  =  (  lb  +  It  )/2  and  1  =  (  lb  -  lt  )/2. 

Fig.  3.23  through  3.32  show  radiation  patterns  for  a  trapezoidal  charge 
distribution  for  different  values  of  top  and  base  lengths.  All  except  the  last 

three.  Fig.  3.30  through  Fig.  3.32  are  for  i ?  =  12.  TFTL  is  the  symbol  for  the 
top  level  length  of  the  trapezoidal  function  and  TFBL,  for  the  base  level 
length  of  trapezoidal  function. 

As  mentioned  earlier,  the  form  factor  is  sensitive  to  the  ratios  of  the 
characteristic  lengths  to  the  wavelength  of  the  emitted  radiation.  Here  the 

appropriate  ratios  are  ^77  and  1/A.  In  Fig.  3.23  both  these  ratios  are  less  than 
one  so  that  the  radiation  pattern  is  characteristic  of  a  point  charge 
distribution.  Fig.  3.24  drawn  for  ?/\  1.2  and  1/A  =  0.29  shows  a  different 
pattern  with  the  form  factor  already  becoming  quite  evident. 

Fig.  3.25  through  Fig.  3.27,  calculated  for  a  constant  slightly  larger  ratio 
of  e/\  =  1.71,  and  increasing  values  of  1/A  are  all  similar  since  the  radiation 
wavelength  is  too  large  to  be  sensitive  to  the  smaller  details  of  the  charge 
distribution.  The  same  comment  applies  also  to  Fig.  3.28  for  which  (/\  is  still 
larger  but  1/A  is  less  than  one. 
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Fig.  3.32  Trapezoidal  Function  with  (/ A  =  1.71.  1/A  =  0.') 


The  last  figure  of  the  q  =  12  series  is  drawn  for  both  (/  X  and  1/A  greater 
than  one.  Consequently  one  would  expect  the  radiation  envelope  to  display 
both  terms  in  the  form  factors  for  a  trapezoidal  charge  distribution.  That  it 
apparently  does  not  may  be  due  to  the  fact  that  the  lengths  chosen  in  order 
to  achieve  these  make  for  a  charge  distribution  which  is  almost  triangular. 

Figures  3.30-3.32  are  calculated  for  the  below  threshold  condition  of  q  = 
37.  All  have  (/X  greater  than  one  and  1/A  less  than  one.  As  with  other  figures 
with  this  value  of  q ,  the  envelope  is  more  evident. 

When  both  values  of  1/ A  and  1/A  increase  the  peak  values  decrease  mainly 
in  the  region  of  angles  between  0°  and  approximately  55°  and  the  radiation 
intensity  is  enhanced  around  90°.  It  is  difficult  to  determine  the  development 
of  the  radiation  pattern  with  respect  to  t/X  and  1/A  because  both  values  can 
change  the  radiation  pattern  at  same  time.  To  analyze  the  radiation  pattern  in 
another  way  the  angle  where  zeros  happen  can  be  investigated  by  calculation. 
Let's  set: 


^  -  n  7T  , 


T  =  m  ’ 


(3.1) 


where  kz  =  k  cos  0,  and  m,n  are  integers.  Substituting  for  ^  gives 


cos  0n  = 


±  n  -y- 


cos  0m  =  ±  m 


(3.2) 


~r 


If  as  in  most  cases  considered  here 


<  1,  these  are  no  solutions  for  0m. 


Thus  only  the  solutions  for  9n  are  considered  below.  A  example  is  Fig.  3.25 

l 

where  —y  =  1.71,  f  =  6  cm,  A  =  3.5  cm. 


cos  9  = 


The  solutions  are  for  9n  are 


9\  =  54.5°  and  9.  i  =  125.5°. 


=  ±  n  *  0.58 


which  agree  with  the  figure.  These  values  also  apply  to  Fig.  3.26.  Fig.  3.2' 

f 

Fig.  3.31,  and  Fig.  3.32  which  have  the  same  ratio  of  — y . 

In  Fig.  3.29  where  — t-  1.29.  1/A  =  1.27, 


cos  9n  = 


cos  9„ 


The  solutions  are  9\  =  39.2°,  9. \  =  129.2°  for  9n  and  9\  =  38°.  9-\  =  128°.  In 
This  case  the  zeros  for  9m  and  9n  are  so  close  together  as  to  be 
indistinguishable  in  the  figure. 
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rv.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  object  of  this  work  is  to  assess  the  effects  of  the  charge  distribution 
of  a  bunch  on  the  sub-Cerenkov  radiation  from  a  charged  particle  beam, 
specifically  investigating  the  short  beam  path  region  at  energies  above  and 
below  the  threshold.  It  contains  many  calculated  radiation  patterns  for  the 
three  kinds  of  charge  distributions. 

With  the  charge  distribution  of  Gaussian  function,  the  radiation  intensity 
envelope  near  90°  represents  the  square  of  the  bunch  form  factor  F(k)  when 
the  bunch  length  is  larger  than  a  half  of  the  wavelength.  When  the  bunch 
length  is  about  same  as  the  wavelength  or  larger  than  that  of  the  level  charge 
distribution  the  radiation  pattern  around  90°  is  also  dominated  by  the  square 
of  the  bunch  form  factor  F(k). 

It  is  difficult  to  determine  the  characteristics  of  the  radiation  pattern  on 
the  trapezoidal  function  by  the  ratio  of  bunch  length  to  wavelength  because 
the  level  top  length  and  level  base  length  can  change  the  shape  of  radiation 
pattern  at  same  time.  However  the  radiation  pattern  of  the  trapezoidal 
function  can  be  partially  characterized  by  analyzing  the  zeros  of  the  curves. 

A  summary  of  results  is  : 

1.  When  the  bunch  length  is  much  smaller  than  the  wavelength  of  the 
emitted  radiation,  the  radiation  pattern  is  characteristic  of  point 
charge. 

2.  When  tj  is  increased  the  total  radiation  pattern  does  not  change 
except  for  minor  differences  in  the  frequency  of  oscillation. 


3. 


When  the  wavelength  of  the  emitted  radiation  is  comparable  to  or 
larger  than  the  bunch  length,  a  significant  amount  of  radiation 
appears  at  near  90°  to  the  beam. 

4.  The  envelope  of  the  90°  radiation  is  the  square  of  the  bunch  form 
factor.  This  envelope  is  skewed  slightly  as  a  result  of  monotonic  character  of 
the  geometrical  envelope  function. 

To  enhance  this  preliminary  work  the  following  is  recommanded: 

1.  Investigate  the  behavior  of  the  square  of  form  factor  F(k)  only  when 
the  bunch  length  varies.  This  will  give  the  insight  into  the  behavior 
of  radiation  pattern. 

2.  Determine  why  varying  the  top  length  of  the  trapezoidal  pulse  shape 
causes  significant  enhancement  of  the  forward  lobe  envelope. 

3.  Determine  the  effects  of  changes  in  beam  energy,  pulse  width. pulse 
frequency,  and  path  length. 

4.  Confirm  the  results  found  by  mapping  the  microwave  radiation 
measured  at  a  comparable  particle  accelerator. 

The  examples  presented  here  are  idealized  and  there  has  been  no  attempt  to 
take  into  consideration  such  factors  as  beam  pulse  erosion  as  it  transverses  the 
atmosphere  and  the  ground  plane  reflection  in  interference[Ref.  4:  p.  1 996] .  But 
the  complete  angular  map  of  the  envelope  of  the  radiation  pattern  could 
provide  a  measure  of  F(kz)  and  an  excellent  reference  to  establish  a  method  of 
determining  the  beam  pulse  charge  distribution  based  on  sub-Cerenkov 
radiation  patterns. 
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APPENDIX  A 


PROGRAM  FOR  CALCULATION  OF  RADIATED  POWER 


*  THIS  PROGRAM  IS  TO  CALCULATE  THE  POWER  PER  UNIT  SOLID  ANGLE  RADIA-  * 

*  TED  AT  THE  FREQUENCY  NU .  IT  WAS  ORIGINALLY  DEVELOPED  BY  PROFESSOR  * 

*  JOHN  R.  NEIGHBOURS  AND  OTHERS  FOR  TECTRONIX  COMPUTER.  THIS  AUTHOR  * 

*  MODIFIED  TO  RUN  ON  IBM  370  MAINFRAME  COMPUTER  IN  DEC.  1987. 

*  THE  THEORY  AND  ALL  EQUATION  ARE  WELL  GIVEN  IN  REF.  5. 


PROGRAM  JNN 

COMMON  / VA/HH , HU , SDXX  /DOUP A/L , ETA , TMAX , DT , E . N , B , LL , RR , NN , RFTL . 
+  RFBL , TFTL , TFBL , TRIBL , NBETA 

DOUBLE  PRECISION  C , C0 ,MU0 , NU0 , Q , P I , E0 , E , NBETA , 

DOUBLE  PRECISION  TMAX . DT , NU , L , ETA . N , B , LL , RR , NN 

DOUBLE  PRECISION  I IMAX , KK , CONST ,WMAX , THETA ,W, T , KKZ , H , DENOM1 . 

DOUBLE  PRECISION  DENOM2.K0 

DOUBLE  PRECISION  GAMMA , V , U . I I , R , HA , HB , HC , HD , DD , DEXIN 
DOUBLE  PRECISION  RFTL, RFBL, TFTL. TFBL. TRIBL 

DOUBLE  PRECISION  DENOM3 , DEN0M4 , DENOM5 , DEN0M6 , DEN0M7 , DEN0M8 , 
DOUBLE  PRECISION  HE.HF.HG.HI 

DOUBLE  PRECISION  HJ , HK , HL , Z , LI , LIN , HX , HW, KE , SDXX 

REAL  XWMAX , XTHETA , XW, YMAX , SDYY , SDY , J  X . J YY , YY , YNUM , X , Y , SDX , 

+  XX  , XNUM , RTMAX , RSDXX , RWMAX , RTHETA , RW 

REAL  REL , REE . REN , REB , RELL . RERR , RENN , DELTA 
REAL  RERFTL , RERFBL , RETFTL , RETFBL , RETRIB , REHU 
INTEGER  I , K , KMAX . IMAX , HH . J , JXMAX , JYMAX , IPA , P , HU 
CHARACTER* 1  P ARACH , EXPPTS , EXPTCH , SAVEQ . PPRE , CHANGE , ANS , AXCH , 

+  ANS2.ANS3 

CHARACTER* 13  NAMEPF , NAMEWF , NAMEEX 
DIMENSION  XTHETA( 1  ; 1 00  ) 

DIMENSION  XW( 1:100) 

DIMENSION  THETA( 0:4500) 

DIMENSION  W( 0:4500) 

DIMENSION  H( 1 : 7 ) 

C 

C  *  ASSIGNMENT  OF  CONSTANTS 

C 

C=2 . 997 1 22D08 

C0  =  2  . 997925D08 

MU0  =  1  .2566437061D-06 

C 

C  . . TO  VARY  CHARGE  ****** . 

r> 

v 

PI =3. 1 41 5926535897932D0 
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’TTTTTTWT’ 


c 

5 


45 


55 

C 


C 

C 

C 


C 


E0=0 . 5 1 17D0 
PRINT*,  • 

PRINT*,  ’  ENTER  BUNCH  CHARGE .( DEFAULT=  1.5E-12)  ’ 

READ*,  Q 

I F ( Q . GT . 0 )  GO  TO  45 
Q= 1 .5D-12 

PRINT*,  ’  ENTER  FUNDAMENTAL  FREQUNCY  NU0.  ( DEFAULT=  2.85E9)' 
READ*,  NU0 

IF( NU0 . GT . 0 )  GO  TO  55 
NU0=2 . 85D09 
CONTINUE 


PRINT*, ’THE  INDEX  OF  REFRACTION  IN  THE  AIR  IS  1  . 000268( DEFAULT )  . 
PRINT*,’ DO  YOU  WISH  TO  CHANGE?  (Y/N)  ’ 

READ  ’ ( 1  A )  ’  ,  ANS2 
IF( ANS2 . EQ . ’ Y’ )  THEN 

PRINT*, ’ENTER  THEINDEX  OF  REFRACTION. 

READ*,  DEXIN 
ELSE 

DEXIN= 1 . 000268D0 
END  IF 
C=C0/DEXIN 


»««**»» 


DECISION  OF  PARAMETER  FILE******* 


PRINT  * 

PRINT  *,  ' 
PRINT  *,  ’ 
PRINT  *,  ’ 
PRINT  * 

PRINT  *, ’ 
READ  ' ( 1  A) '  , 
IF( PARACH . EQ 


1:  CREATE  NEW  PARAMETER  FILE’ 
2:  USE  OLD  FILE’ 

3:  EXIT’ 


ENTER  YOUR  CHOICE 
PARACH 

'3' )  GOTO  700 
IF( PARACH. EQ. ’2' )  CALL  INAMPF 
IF( PARACH. EQ. ’ 1 ’ )  CALL  MAKEPF 


READ  IN  ALL  PARAMETERS  FROM  EXTERNAL 
FILE/PRINTOUT  PARAMETERS 


C 


99 

WRITE( 6,100) 

1  00 

FORMAT(/// ’THE  PARAMETERS  YOU  ENTERED 
REWIND  1 

ARE:  '  ) 

120 

READ( 1 ,*)  HH 

131 

WRITE( 6,135)  HH 

1  35 

FORMAT( / ’FOURIER  TRANSFORM  CHOICE  : 

IF  (HH.NE. 1  )  GOTO  1  37 

WRITE  (6,136) 

’  ,18) 

136 

FORMAT( / ,37X, ’GAUSSIAN  FUNCTION  ’  ) 

137 

IF  (HH.NE. 2)  GOTO  139 
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V  > 


V  ». 


WRITE  (-=■ ,  138) 

FORMAT  (/,37X.  ’LEVEL  FUNCTION  ’  ) 

IF  (HH.NE.3)  GOTO  141 
WRITE  (6,140) 

FORMAT  (/’SUM  OF  LEVEL  AND  RIPPLE  FUNCTIONS'  ) 
IF  (HH.NE.4)  GOTO  143 
WRITE  (6,142) 

FORMAT(/ ’MULTIPLE  HUMP  FUNCTION'  ) 

IF  (HH.NE.5)  GOTO  145 
WRITE  (6,144) 

FORMAT  (/’ROUNDED  FUNCTION'  ) 

IF  (HH.NE.6)  GOTO  147 
WRITE  (6,146) 

FORMAT  (/’TRAPEZOIDAL  FUNCTION’  ) 

IF  (HH.NE.7)  GOTO  149 
WRITE  (6,148) 

FORMAT  (/’TRIANGLE  FUNCTION’  ) 


,  D1  6 . 4  ) 


READ( 1 , * )  ETA 
WRITE( 6,151 )  ETA 

FORMAT ( / ’  ETA  ( L/RAMMDA )  :  ’.D16.4) 

READ( 1 , * )  TMAX 
WRITE( 6,155)  TMAX 

FORMAT ( / '  MAXIMUM  THETA  VALUE  :  ’.D16.4) 

READ( 1 , * )  DT 
WRITE( 6,165)  DT 

FORMAT ( / ’  THETA  CALCULATION  INCREMENT  :’,D16.4) 
READ( 1 ,*)  SDXX 
WRITE( 6,175)  SDXX 

FORMAT ( / '  THETA  MARKING  INCREMENT  :  ’,D16.4) 

READ( 1 ,*)  E 
READ( 1 , * )  NBETA 
WRITE( 6,185)  E 
WRITE( 6,185)  NBETA 

FORMAT ( / '  BEAM  ENERGY  :  ’.D16.4 

FORMAT ( / '  VALUE  OF  NBETA  :  ’,D16.4 

READ( 1 ,*)  N 
WRITE( 6,195)  N 

FORMAT ( / '  HARMONIC  NUMBER  :  ’,D16.4) 


,  D1 6 . 4  ) 


’  ,  D1  6 . 4  ) 
’  ,  D1  6 . 4  ) 


,D1 6.4 


IF  (HH.NE. 1 )  GOTO  206 


READ( 1 , * )  B 
WRITE( 6,205)  B 

FORMAT ( / ’  BUNCH  LENGTH  (M)  : 


’  ,  D1  6 . 4  ) 


IF  (HH.LT.2)  GOTO  216 
IF  (HH.GT.4)  GOTO  216 


READ( 1 , * )  LL 
WRITE( 6,215)  LL 


.  ^  HJ f  I,**.**,  ^ ■^‘«1  ■*,  •*,  »*,  »*«  ■’  .  •*  »*  •*,  *'  **  *" 


215  FORMAT ( / '  LEVEL  FUNCTION  LENGTH  (M)  :  ’,316.4) 

C 

216  IF  (HH.NE.3)  GOTO  9235 
C 

READ( 1 , * )  RR 

224  WRITE( 6 , 225 )  RR 

225  FORMAT ( / ’  FOURIER  RATIO  CONSTANT  :  ’.D16.4) 

C 
C 

READ( 1  ,*)  NN 

234  VRITE( 6 , 235 )  NN 

235  FORMAT( / '  NUMBER  OF  RIPPLES  :  ’,316.4) 

C 

9235  IF  (HH.NE.4)  GOTO  236 
C 

READ  (1 , *)HU 

9237  WRITE  (6,9239)  HU 

9239  FORMAT( / 1  NUMBER  OF  HUMPS  :  '.14) 

C 

236  IF  (HH.NE.5)  GOTO  241 
C 

READ  ( 1 , * )  RFTL 

9238  WRITE  (6,238)  RFTL 

238  FORMAT  (/’  ROUNDED  FUNCTION  LEVEL  TOP  LENGTH  :  ’.D16.4) 

C 

READ  ( 1 , * )  RFBL 

9240  WRITE  (6,240)  RFBL 

240  FORMAT  (/>  ROUNDED  FUNCTION  BASE  LENGTH  :  ’.D16.4) 

C 

241  IF  (HH.NE.6)  GOTO  246 
C 

READ  ( 1 , * )  TFTL 
9243  WRITE  ( 6 , 243  )TFTL 

243  FORMAT  (/’  TRAPEZOIDAL  FUNCTION  TOP  LEVEL  LENGTH  :’,D16.4) 
C 

READ  ( 1 , * )  TFBL 
9245  WRITE  (6,245)TFBL 

245  FORMAT  (/’  TRAPEZOIDAL  FUNCTION  BASE  LENGTH  :  ’ .D16.4) 

C 

246  IF  (HH.NE.7)  GOTO  249 
C 

READ  ( 1 , * )  TRIBL 
9248  WRITE  (6,248)TRIBL 

248  FORMAT  (/’  TRIANGLE  FUNCTION  BASE  LENGTH  :  ’.D16.4) 

249  WRITE( 6,250) 

250  FORMAT( / / ’ DO  YOU  WISH  TO  CHANGE  ANY  OF  THESE  ?(Y  OR  N)') 
READ’ ( A1 ) ' ,  CHANGE 

C 

IF( CHANGE. EQ. ’ Y ’ )  GO  TO  700 
C 

*»«**»«»«»##«««»#««»»»»«»»«»»«» 

*  ***  BEGIN  CALCULATION  ***  * 


I IMAX=TMAX/DT 
IMAX=INT(  UMAX) 

NU=N*NU0 

IF  ( ( HH . GE .  2 ) . AND . ( HH . LE . 4 ) )  THEN 
DD=LL/2 . 0D0 
END  IF 

KK=2.0D0*PI*NU/C 

CONST =MU0*C*( NU0**2 .0)*(Q**2.0)/(8. 0D0#( PI**2 . 0 ) ) 
WPIAX  =  0 . 0D0 
DO  750  1=0, IMAX 

THETA( I )=DBLE( I ) *DT 
T=PI*THETA( I )/ 1 80 . 0D0 
KKZ=KK*DCOS(T) 

IF  (HH.NE. 1 )  GOTO  326 

H( 1 )=DEXP( -1 . 0D0*( ( KKZ*B/2 . 0D0 )**2 . 0 ) ) 

CONTINUE 

IF  (HH.LT.2)  GOTO  328 
IF  (HH.GT.3)  GOTO  328 

DENOM1 =KKZ*DD 

IF  ( DENOM1 . EQ . 0 )  GOTO  327 

H( 2  )  =  ( DSIN  (KKZ*DD) )/(KKZ»DD) 

GOTO  328 
H( 2 ) =0 . 0D0 

IF( HH . NE . 3 )  GOTO  333 

KO=PI*NN/DD 
DEN0M2= ( KKZ-KO )*DD 

IF( DEN0M2 . EQ . 0 )  GOTO  329 

HA=  DSIN( DEN0M2 ) /DENOM2 

GOTO  330 

HA=0 . 0D0 

DENOM3= ( KKZ+KO ) *DD 
IF  ( DENOM3 . EQ . 0 )  GOTO  331 
HB=DSIN( DENOM3 )/DENOM3 
GOTO  332 
HB=0 . 0D0 


H(3)=H(2)+((RR/2.3D0)*( HA+HB )  ) 


IF( HH . NE . 4  )  GOTO  335 

Z=(KKZ*DD) 

IF( Z . EQ . 0 )  GOTO  334 

H( 4 ) = 1 . 0D0* ( (DSIN(Z) )/Z) 

DEN0M4=  1 .0 
DO  9333  P  =  1 , HU 

KE= ( 2 . 0D0*DBLE( P ) ) - 1 .0D0 

LI  =  ( ( 2 . 0D0*DBLE( HU ) ) -KE ) / ( ( 2 . 0D0*DBLE( HU ) ) -  1  . 0D0  ) 

LIN=( ( -1  . 0D0 )**DBLE(P ) )*( -1 . 0D0  ) 

HX=(LIN*LI*DSIN(Z*LI ) )/(Z*LI ) 

HW=LIN*LI 
H(4)=H(4)+HX 
DEN0M4=  DEN0M4+HW 
CONTINUE 
H(4)=H(4 ) /DEN0M4 

GOTO  335 

H( 4 ) =0 . 0D0 

IF  (HH.NE.5)  GO  TO  339 

HL=( 8D0* 1 6D0 ) / ( ( ( RFBL-RFTL ) *  *2 . 0 ) * ( KKZ * *3 ) * ( RFBL+RFTL  )  ) 
HG=HL*DSIN( ( KKZ/ 4 . 0D0 ) *( RFEL+RFTL ) ) 

H J  =  ( RFTL*DSIN( KKZ*RFTL/2 . 0D0 ) ) / ( KKZ*RFTL/2 . 0D0 ) 

HK=  ( RFBL*DSIN( KKZ-RFBL/2 . 0D0  ) ) / ( KKZ*RFBL/2 . 0D0  ) 

HI  =  ( 32D0 / ( ( ( ( RFBL-RFTL ) * ( KKZ ) ) *  *2 . 0 ) * ( RFTL+RFBL ) ) ) * ( HJ+HK ) 
H( 5 ) =HG-HI 

IF  (HH.NE.6)  GOTO  345 

DENOM6=KKZ* ( TFBL+TFTL ) /4 . 0D0 
DEN0M7 =KKZ* ( TFBL-TFTL ) /4 . 0D0 
HE=DSIN( DENOM6 ) /DENOM6 
HF=DSIN( DEN0M7 ) /DEN0M7 
H( 6  )  =  ( TFBL+TFTL ) /2 . 0D0*HE*HF 

IF( HH . NE . 7 )  GOTO  349 

DENOM5=(KKZ*TRIBL)/4.0D0 

IF  (DENOM5.EQ.0)  GOTO  346 
H( 7 )  =  ( DSIN( DENOM5 ) / ( DENOM5 ) ) **2 . 0 

GOTO  349 


H( 7 ) =0 . 0D0 


CONTINUE 


v.  v 


GAMMA =E/E0 
L=ETA*C/NU 

GAMMA  =  DSQRT( 1 .000/(1 . 0D0- ( NBETA/DEXIN )  *#2 ) ) 

V  =  C0*(DSQRT(  1  . 0D0-(  1  . 0D0/ ( GAMMA**2 . 0 ) ) ) ) 

V=C0*( NBETA/DEXIN ) 

U=(KK*L/2.0D0)*( (C/V)-DCOS(T)) 

1 1 =DSIN( U  )/U 
R=KK*L*H(HH)*II*DSIN(T) 

W( I ) = CONST * ( R**2 . 0 ) 

CALCULATION  OF  MAXIMUM  V  *** 

WMAX=DMAX1 ( WMAX , W( I ) ) 

CONTINUE 

***  DECISION  TO  SAVE  VALUES  CALCULATED  . * 

WRITE( 6 , 360  ) 

FORMAT ( / ’ DO  WANT  TO  SAVE  THE  RESULTS  IN  A  FILE?  (Y/N)’  ) 
READ1 (A1 ) ' ,  SAVEQ 

IF(SAVEQ.EQ.  ’N’ )  GOTO  4-40 

PRINT  *.  'UNDER  WHAT  FILE  NAME?  ’ 

READ  1 ( A 1 0 ) ' ,  NAMEWF 

OPEN(UNIT=3 , FILE=NAMEWF , STATUS* ' NEW  -  , FORM= ’ FORMATTED '  ) 

DO  370  1=0, IMAX 

RTHETA  =  REAL( THETA( I )  ) 

RW=REAL( W( I ) ) 

WRITE( 3 . 380 )  RTHETA.RW 
FORMAT(4X,F8.3,2X,F20. 18) 

CONTINUE 
WRITE( 3 , 385  ) 

FORMAT( 7X , '0.000' ) 

WRITE( 3 , 390 )  WMAX 
F0RMAT(2X. 'WMAX  =  '.D16.4) 

WRITE( 3 , 392 )  ETA 

FORMAT ( 2X , ’ ETA  =  ’.D16.4) 

WRITE( 3 , 394 )  NBETA 
FORMAT( 2X , ' NBETA  =  ’.D16.4) 

WRITE( 3 , 396  )  B 

FORMAT(2X, ’ BUNCH  LENGTH  =  '.D16.4) 


***  DECISION  TO  PRINTOUT , RE-RUN , OR  EXIT 
WRITE( 6 , 445 ) 

FORMAT(////’  1:  PRINT  OUT  VALUES’,/ 

’  2:  RUN  PROGRAM  AGAIN’ ,/ 


w  rTj 


3:  EXIT’,/// 

-  ’  ENTER  CHOICE. _ ’ ) 

READ ' ( A1  ,  '  ,  PPRE 

IF' PPRE . EQ.  ’ i  ’  )  GOTO  66a 
IF' PPRE. EQ. ’2  ’  )  GOTO  05 
IF1 PPRE . EQ . ’3’ )  GOTO  700 
GOTO  440 

C 

C  •**  PRINT  OUT  VALUES  *»* 


660  WRITE! 6 , 670 ) 

673  FORMAT!’  THETA  W/23(  ’-*)/) 

DO  690  1=0, IMAX 

WRITE! 6, 680)  THETA! I ) ,W( I ) 

680  FORMAT! F8. 2, FI  5. 4) 

690  CONTINUE 
C 

. . *  PROGRAM  ENDS . ****.„*,*,*.**«* 

700  END 

C 
C 
C 


*  ***  SUBROUTINE  INAMPF  ***  * 

*«***»«**»«»*»#**»***«#«***»» 

C 

SUBROUTINE  INAMPF 

C 

C  ***  INPUT  OLD  PARAMETER  FILE  NAME 
CHARACTER* 10  NAMEPF 
WRITE! 6 , 900 ) 

900  FORMAT!/ ’ENTER  THE  NAME  OF  THE  PARAMETER  FILE  YOU  WISH  TO’ 
+  ’  USE! MAX  10  CHAR)  ’  ) 

READ ’ ( A1 0 ) * ,  NAMEPF 
WRITE! 6,91 1  ) 

911  FORMAT! / ’TEST  LINE  ONE’  ) 

OPEN! UNIT= 1 , FILE=NAMEPF , STATUS* ’ OLD ’ , FORM* ’ FORMATTED ’ , 

+  BLANK* ’ NULL ’  ) 

WRITE! 6,912) 

912  FORMAT! / ’TEST  LINE  TWO’  ) 

RETURN 

END 

C 


»  **»  SUBROUTINE  INAMEX  *»*  * 


C 

SUBROUTINE  INAMEX 
C 

C  ***  INPUT  OLD  EXPERIMENTAL  POINTS  FILE  NAME  »** 
CHARACTER* 10  NAMEEX 


WRITE(6, 1000) 

1000  FORMAT( / 'ENTER  NAME  OF  EXPERIMENTAL  POINTS  FILE  YOU  WISH’, 

+  -TO  USE( MAX  10  CHAR)  '  ) 

READ'  ( A1  0 )  ’  ,  NAMEEX 

OPEN( UNIT  =  2 , FILE  =  NAMEEX , STATUS = ’ OLD 1 ,FORM= 'FORMATTED' , 

+  BLANK= ’ NULL ’  ) 

RETURN 

END 

C 

C 

*###*#**##***#**#■***■»******** 

*  ***  SUBROUTINE  MAKEPF  ***  * 

#*##*#*#*#*»****»********»»*# 

C 

•  SUBROUTINE  MAKEPF 
C 

C  ***  INPUT  OF  PARAMETERS  INTO  NEW  FILE  (UNIT=1)*** 

C 

COMMON  /VA/HH , HU , SDXX  /DOUPA/L , ETA , TMAX , DT , E , N , B , LL , RR , NN , RFTL , 
+  RFBL , TFTL , TFBL , TRIBL , NBETA 

DOUBLE  PRECISION  L , ETA , TMAX , DT , E , N , B , LL , SDXX . NBETA , 

+  RR , NN , RFTL , RFBL . TFTL , TFBL .TRIBL 

INTEGER  HH.HU 
CHARACTER* 10  NAMEPF 
WRITE( 6 , 1 300 ) 

1300  FORMAT(/ ’ENTER  THE  NAME  FOR  THE  NEW  PARAMETER  FILE.’, 

+  ’(MAX  10  CHAR)  ’  ) 

READ ’ ( A1 0 ) ’ ,  NAMEPF 

OPEN( UNIT= 1 , FILE=NAMEPF, STATUS*’ NEW ’, FORM* ’FORMATTED' , 

+  ACCESS- ' SEQUENTIAL ’ . BLANK* ’ NULL ’ ) 

PRINT*.  ’ 

WRITE( 6 , 1 320 ) 

1320  FORMAT(///// ’NOTE:  ALL  THE  PARAMETERS  MUST  BE  ENTERED', 

+  ’  ALTHOUGH  SOME  MAY  NOT  BE  USED’) 

WRITE( 6 , 1 330 ) 

1330  FORMAT ( / '  1:  GAUSSIAN  FUNCTION’/’  2:  LEVEL  FUNCTION’, 


+ 

/• 

3: 

SUM  OF  LEVEL  AND  RIPPLE  FUNCTIONS 

+ 

/’ 

4: 

MULTIPLE  HUMP  FUNCTION’, 

+ 

r 

5: 

ROUNDED  FUNCTION’, 

+ 

r 

6: 

TRAPEZOID  FUNCTION’, 

+ 

/' 

7: 

TRIANGLE  FUNCTION’ ) 

WRITE( 6 , 1  340  ) 

1340  FORMAT( / 'ENTER  THE  FOURIERTRANSFORM  CHOICE:’) 

READ*,  HH 
C 

1356  I F ( HH . LE . 7 )  GO  TO  1359 
C 

WRITE( 6 , 1 357 ) 

1357  FORMAT ( ’  FOURIER  TRANSFORM  CHOICE  NOT  DEFINED.’///, 

+  ’  ENTER  NEW  CHOICE: ’ ) 

READ*,  HH 
C 


c 

1359  WRITE( 1 ,*)  HH 
VRITE(6, 1 370) 

1370  FORMAT( / 'ENTER  ETA  (L/RAMMDA).  (  DEFAULT  =  12.0D-2  )  '  ) 

READ*,  ETA 

IF( ETA . GT . 0 . 0 )  GOTO  1385 
ETA= 1 2 . 0D0 
1385  VRITEO  .*)  ETA 
¥RITE( 6 , 1 400  ) 

1400  FORMAT( / 'ENTER  THE  MAXIMUM  VALUE  OF  THETA  TO  BE', 

+  '  CALCULATED(DEG) . ' / '  (  DEFAULT  =  45-0  ’) 

READ*,  TMAX 

IF( TMAX . GT . 0 . 0 )  GOTO  1415 
TMAX=45 . 0D0 
1415  WRITE( 1 ,*)  TMAX 
WRITE(6, 1430) 

1430  FORMAT( / ’ENTER  THE  ANGLE  CALCULATION  INCREMENT  DEG) ’  , 

+  ' ( DEFAULT  =  0 . 1  )  '  ) 

READ*.  DT 

IF( DT . GT . 0 . 0 )  GOTO  1445 
DT  =  0 . 1 0D0 

1445  WRITE( 1 ,*)  DT 
WRITE(6, 1460) 

1460  FORMAT( / 'ENTER  THE  ANGLE  MARKING  INCREMENT  DEG ) ’  , 

+  ’  ( DEFAULT  =  1 5  ■  0  )  '  ) 

READ*,  SDXX 

IF( SDXX . GT . 0 . 0 )  GOTO  1475 
SDXX= 1 5  - 0D0 
1475  WRITE(1 ,*)  SDXX 
WRITE( 6 , 1 490 ) 

Cl  490  FORMAT(/ ’ENTER  THE  BEAM  ENERGY  ( MEV )  ( DEFAULT= 100.0)  ’  ) 
1490  FORMAT(/ 'ENTER  THE  VALUE  OF  NBETA  ( DEFAULT  =  0 . 98 )  ’  ) 

C  READ* ,  E 

READ*,  NBETA 

C  I F ( E . GT . 0 . 0 )  GOTO  1505 

IF( NBETA. GT. 0.0)  GOTO  1505 
C  E= 1 00 . 0D0 

E=0 . 98D0 

Cl  505  WRITE( 1  ,*)  E 
1505  WRITE( 1 ,*)  NBETA 
WRITE( 6 , 1 520 ) 

1520  FORMAT( / 'ENTER  THE  HARMONIC  NUMBER  ( DEFAULT= 1.0)  ’  ) 

READ*,  N 

IF( N . GT . 0 . 0 )  GOTO  1535 
N= 1 . 0D0 

1535  ¥RITE( 1 ,*)  N 
C 

IF  (HH.NE. 1 )  GOTO  1575 
C 

WRITE( 6,1550) 

1550  FORMAT( / 'ENTER  THE  BUNCH  LENGTH  (M).(  DEFAULT  =  0.24D-2  ) 


READ*.  3 
C 

IF( B . GT . 0 . 0 )  GOTO  1565 
C 

B=0 . 24D0 

1565  WRITE( 1 ,»)  B 
C 

1575  IF  ( ( HH . LT . 2 ) . OR . ( HH . GT . 4 ) )  GOTO  1605 
C 

VRITE( 6 , 1 580  ) 

1580  FORMAT(/ 'ENTER  THE  LEVEL  FUNCTION  LENGTH  (M).’. 

+  '  ( DEFAULT  =  0 . 01  )  '  ) 

READ*.  LL 
C 

IF( LL . GT . 0 . 0 )  GOTO  1595 
C 

LL=0 . 01 D0 

. ..5  WRITE( 1 ,*)  LL 

C 

1605  IF  (HH.NE.3)  GOTO  1635 

C 

WRITE( 6,1610) 

1610  FORMAT(/ ///'ENTER  THE  FOURIER  RATIO  CONSTANT  ( DEFAULT* 1.0).') 

READ*,  RR 
C 

IF( RR . GT . 0 . 0 )  GOTO  1625 
C 

RR= 1 . 0D0 

1625  WRITE( 1 ,*)  RR 
C 

1635  IF  (HH.NE.3)  GOTO  1665 
C 

WRITE(6, 1640) 

1640  FORMAT( / 'ENTER  THE  NUMBER  OF  RIPPLES  ( DEFAULT= 1.0)  ’  ) 

READ*.  NN 
C 

IF( NN . GT . 0 . 0 )  GOTO  1655 
C 

NN= 1 . 0000 

1 655  WRITE( 1 , * )  NN 
C 

1665  IF  (HH.NE.4)  GOTO  1700 
C 

WRITE( 6,1670) 

1670  FORMAT(  / 'ENTER  NUMBER  OF  HUMPS  (DEFAULTS)  '  ) 

READ*,  HU 
C 

IF( HU . GT . 2 )  GOTO  1680 
C 

HU  =  2 

1680  WRITE(  1  ,*)HU 
C 
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1700  IF  (HH.NE.5)  GOTO  1800 
C 

WRITE  (6,1710) 

1710  FORMAT{/ ’ENTER  THE  ROUNDED  FUNCTION  LEVEL  TOP  LENGTH  (CM)', 
+  / '  (DEFAULT  =  1.0)  ’  ) 

READ*,  RFTL 
C 

IF  ( RFTL. GT. 0.0)  GOTO  1730 
C 

RFTL  =  1 . 0D0 
1730  WRITE  (1,*)  RFTL 
WRITE  (6,1760) 

1760  FORMAT  (/’ENTER  THE  ROUNDED  FUNCTION  BASE  LENGTH  (CM).’/, 

+  ’  (DEFAULT  =2.0)  *  ) 

READ*,  RFBL 
C 

IF  ( RFBL. GT. 0.0)  GOTO  1780 
C 

RFBL  =  1 . 0D0 
1780  WRITE  ( 1 ,*)  RFBL 
C 

1800  IF  (HH.NE.6)  GOTO  1900 
C 

WRITE  (6,1810) 

1810  FORMAT  (/’ENTER  THE  TRAPEZOIDAL  FUNCTION  LEVEL  TOP’, 

+  'LENGTH  (CM).’/’  (DEFAULT  =  1.0)'  ) 

READ*.  TFTL 
C 

IF  ( TFTL. GT. 0.0)  GOTO  1830 
C 

TFTL  =  1 . 0D0 
1830  WRITE  ( 1 ,*)  TFTL 
WRITE  (6,1860) 

1860  FORMAT  (/’ENTER  THE  TRAPEZOIDAL  FUNCTION  BASE  LENGTH  (CM)’, 
+  / ’  (DEFAULT  =2.0)  ’  ) 

READ*,  TFBL 
C 

IF  ( TFBL. GT. 0.0)  GOTO  1880 
C 

TFBL  =1.0 

1880  WRITE  ( 1 ,* )  TFBL 
C 

1900  IF  (HH.NE.7)  GOTO  1990 
C 

WRITE  (6,1910) 

1910  FORMAT  (/’ENTER  TRAI ANGLE  FUNCTION  BASE  LENGTH  (CM).'//, 

+  ’  (DEFAULT  =  1.0)  ’  ) 

READ*,  TRIBL 
C 

IF  ( TRIBL. GT. 0.0)  GOTO  1930 
C 

TRIBL  =  1 . 0D0 
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*  * 


*  THIS  PROGRAM  IS  TO  PLOT  THE  POWER  PER  UNIT  SOLID  ANGLE  RADIATED 

*  AT  THE  FREQUENCY  NU.  THE  THEORY  AND  ALL  EQUATIONS  ARE  DESCRIBED 

*  IN  REF. 3.  IT  ASSUMES  THETA  IN  COL.  1  (DEGREES)  AND  POWER  IN  COL.  2  * 

*  (WATT)  RESPECTIVELY.  IT  SUPPORTS  A  SINGLE  PLOT  OF  MULTIPLE  CURVES  * 

*  AND  USES  THE  COMPRS/DISSPOP  OPTION  OF  DISSPLA. 


C 

REAL  YX , XMIN , XMAX , YMIN , YMAX , X , Y , DX , AX , AY , XPHY , YPHY 
INTEGER  I OUT , ITICK, ITICKX , ITICKY , I YVAR , ICT , INN , NB ,NETA , NF , 

+  NBL.LF.NMOD, NRAM.TLL.BL.NL 

CHARACTER* 1  YES, ANS 
CHARACTER**  BLNK( 2 ) 

CHARACTER*60  IMESS , IPHI , I RAM , IMOD , INB , IBL , I ETA , LFL . INF , ITLL , 
+  INL 

C 

DIMENSION  X( 2000  )  , Y( 2000 ) 

DIMENSION  YX( 1 : 5000 , 1 : 2 ) 

C 

DATA  YES/’Y'/ 

DATA  BLNK/ 1  ' , ’  ' / 

EPS= .001 
C 

C  INITIALISE  PLOTTER 

C 

100  CONTINUE 

9400  FORMAT ( A1 ) 

9401  FORMAT ( I 1 ) 

9402  FORMAT ( 12  ) 

WRITE( 4 , 9403  ) 

9403  FORMAT ( ’  DO  YOU  WANT:',/, 

X  ’  1 .  COMPRES’ ,/ , 

X  ’  2.  TEK6 1 8 ' ,/. 

X  ’  ENTER  NUMBER  (II)’,/) 

READ( 4 , * )  IOUT 
IF( IOUT.EQ. 1 )  THEN 
CALL  COMPRS 
ELSEIF( IOUT.EQ. 2 )  THEN 
CALL  TEK618 
ENDIF 

AX  =  4.6810 
AY  =  6.7416 


CALL  HWROT( ’MOVIE’ ) 

CALL  PAGE(8.5, 1 1.0) 

XPHY  =  2. 084-0 
YPHY=2 . 7953 

-/ 

2  DEFINE:  LOWER/UPPER  CASE  STD  AND  GREEK  SHIFT  LABELS  (+,!,?,*) 
2  SHIFT  FOR  INSTRUCTION  SET(A).  AND  THE  TYPE  FACE. 

"i 

k/ 

CALL  MX1 ALF( ’STANDARD’ . ’ • ’ ) 

CALL  MX2ALF( ’L/CSTD’ , ’+’ ) 

CALL  MX3ALF( ’ GREEK 

CALL  MX4ALF( ' L/ CGREEK ’ . ) 

CALL  MX5ALF( ’INSTRU’ , ’&’ ) 

CALL  SWISSL 

CALL  SHDCHR( 90 . . 1 . 0 • 002 , 1  ) 

9410  FORMAT( D1 6.8) 

CALL  HEIGHT ( 0.14) 

CALL  PHYSOR(XPH"  "  HY) 

CALL  AREA2D( AX , AY ) 

C 

C  CLOSE  GRAPH  AND  REOPEN  FOR  LABEL  AXES  NUMBER  TYPES 

C 

CALL  ENDGR( 0 ) 

CALL  PHYSOR(XPHY.YPHY) 

CALL  AREA2D( AX, AY) 

CALL  XNONUM 

C 

C  DETERMINE  THE  AXIS  LABELS: 

C 

C  X-AXIS  LABEL  AND  NUMBER  TYPE 

C 

CALL  EXCMS( 'CLRSCRN' ) 

WRITE( 4,9415) 

9415  FORMAT ( ’  DO  YOU  WANT  INTEGERS  ON  THE  X-AXIS?  (Y/N )’,/) 

READ( 4 , 9400  )  ANS 

IF( ANS.EQ. YES)  CALL  XINTAX 

CALL  XNAME( ’ A+NGLE!  (?Q!,  +DEGREES !  )$  ’ ,100) 

C 

C  Y-AXIS  LABEL  AND  NUMBER  TYPE 

C 

WRITE( 4,9417) 

9417  FORMAT ( ’  CHOOSE  THE  Y  VARIABLE  TO  BE  PLOTTED  (II):'./. 

X  ’  1  -  W  (WATT)  ’  .  / . 

X  *  2 -  DJ  ’  ,  /  ) 

READ( 4 , * )  IYVAR 
IF  (IYVAR.EQ.2)  THEN 

CALL  YNAME( ’  D+J  !$’,100) 

ELSE 

CALL  YNAME('W  ( W+ATT ! ) $ ’ , 1 00 ) 

END  IF 
C 

C  ROTATE  NUMBERS  ON  Y-AXIS  AND  REVERSE  TICKS 


I 
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CALL  YAXANG(0.) 

CALL  XREVTK 
CALL  YREVTK 
VRITE(4,7900) 

7900  FORMAT ( ’  INPUT  YMIN  OF  THE  FIRST  SUBPLOT  (FI  6.8):’,/) 
READ( 4 , *  )  YMIN 
CALL  YAXEND  (  *  NOLAST ’  ) 

YMAX=YMIN+. 1 

CALL  GRAF( 0 . . 1 . , 1 . , YMIN , YMAX , . 1 ) 

CALL  THKFRM( . 027 ) 

CALL  FRAME 
CALL  ENDGR( 0 ) 

C 

C  CLOSED  GRAF  FOR  AXIS  LABELS,  NOW  REOPEN  TO  PLOT 

C 

CALL  PHYSOR( XPHY , YPHY ) 

CALL  AREA2D( AX , AY ) 

CALL  RESET  ( ’ YAXEND ’  ) 

CALL  RESET  (’XNONUM') 

CALL  YNAME( BLNK , 2 ) 

CALL  XNAME( BLNK , 2 ) 

CALL  ENDGR( 0 ) 

C 

C  POSITION  OR  REPOSITION  ORIGIN  AND  SCALE  PLOT 

C 

CALL  PHYSOR( XPHY. YPHY) 

CALL  AREA2D( AX , AY ) 

CALL  RESET  (’YAXEND') 

CALL  YAXEND( ’NOFIRST' ) 

C 

C  DEFINE  COORDINATE  SYSTEM 

C  KEEP  OR  CHANGE  CONSTANTS 

C 

CALL  EXCMS( ’CLRSCRN’ ) 

WRITE( 4 , 9423 ) 

9423  FORMAT ( ’  INPUT  XMIN:  (FI  6. 8)’) 

READ( 4 , * )  XMIN 

WRITE( 4 , 9424 ) 

9424  FORMAT (  '  INPUT  XSTEP :  (FI  6.8)’) 

READ( 4 , * )  XSTEP 

WRITE( 4 , 9425 ) 

9425  FORMAT( '  INPUT  XMAX :  (FI  6.8)’) 

READ( 4 , * )  XMAX 

WRITE( 4 ,9427 ) 

9427  FORMAT ( ’  INPUT  YMIN:  (F16.8)’) 

READ( 4 , * )  YMIN 
WRITE( 4 , 9426 ) 

9426  FORMAT (  '  INPUT  YSTEP :  (FI  6.8)') 

READ( 4 , * )  YSTEP 

WRITE( 4 , 9429 ) 

9429  FORMAT ( ’  INPUT  YMAX:  (FI  6.8)’) 


'.^7^1 


READ(  4- ,  *  )  YMAX 
CALL  EXCMS( ’ CLRSCRN 1 ) 

WRITE(  4- ,  9330  ) 

9330  FORMAT ( ’  HOW  MANY  TICKS  PER  X-AXIS  STEP?  (12)’) 
READ( 4 ,  * )  ITICKX 

CALL  XTICKS( ITICKX) 

WRITE( 4 , 9331 ) 

9331  FORMAT ( 1  HOW  MANY  TICKS  PER  Y-AXIS  STEP?  (12)') 
READ( 4 ,  * )  ITICKY 

CALL  YTICKS( ITICKY) 

IF( XSTEP . EQ . 0 . 0 )  XSTEP = ( XMAX-XMIN ) / AX 
IF( YSTEP . EQ . 0 . 0 )  YSTEP = ( YMAX-YMIN ) / AY 
CALL  GRAF( XMIN , XSTEP , XMAX , YMIN , YSTEP , YMAX ) 

CALL  YNONUM 

CALL  YGRAXS( YMIN, YSTEP, YMAX, AY,  ’  $ ’  , - 1 00 , AX , 0 . 0  ) 
CALL  RESET  (’YNONUM’) 

CALL  XNONUM 

CALI  Y~RAXS( XMIN, XSTEP, XMAX, AX, ’  S’. -100, 0.0, AY) 
CALL  RESET  ( ’ XNONUM ’ ) 

C 

C  CURVE  PLOTTING 

C 

CALL  LINESP ( 1 .8) 

CALL  SPLINE 
CALL  THKCRV( 0.02) 

CALL  GRACE  (0.00) 

C 

C  DRAW  A  CONNECTED  CURVE 

C 

C  MULTIPLE  CURVE  RETURN  POINT 
2000  CONTINUE 

CALL  EXCMS( ’CLRSCRN’ ) 


WRITE( 4 

,9445) 

FORMAT ( 

'  YOU  HAVE  A  CHOICE 

OF  LINE  TYPE: ’ , / , 

+ 

•  1  - 

SOLID  LINE’ , / , 

■f 

’  2  --  -  -- 

DOTTED  LINE’ , / , 

+ 

•  }  - 

DASHED  LINE’ , / , 

+ 

’  4  - 

CHAIN  DOT  LINE’ , / ,  ! 

+ 

'  5  -- - 

CHAIN  DASH  LINE’ , / ,  ' 

+ 

’  PLEASE  INPUT  LINE 

TYPE  NUMBER:  (  11  ) '  ,/ )  ! 

READ(  4  ,  * )  ICT 
IF( ICT . EQ . 2 )  CALL  DOT 
IF( ICT . EQ . 3 )  CALL  DASH 
IF( ICT . EQ . 4  )  CALL  CHNDOT 
IF( ICT . EQ . 5  )  CALL  CHNDSH 
CALL  EXCMS(  ’CLRSCRN’  ) 

C 

C  DATA  INPUT  FROM  EXTERNAL  DATA  FILE 

C 

WRITE(  4- ,  9455  ) 

9455  FORMAT ( ’  INPUT  THE  FILE  DATA  SET  NUMBER  (12)’) 
READ(4,*)  INN 


$ 


80 


wT?j£***j*  v  * 
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NUCOLS  =  2 
IN=0 

DO  2400  11=1 ,5000 

READ( INN, 9460)  ( YX( 1 1 , J J  ) , J J  =  1  .NUCOLS) 

9460  FORMAT ( 4X , F8 . 3 , 2X , F20 . 18) 

IF( ( YX( 11,1)  .LE.  0.0). AND . ( I I . GT . 2 ) )  GO  TO  8000 
IN= IN+1 

2400  CONTINUE 
8000  CONTINUE 
ICOL  =  2 
DO  8010  11=1 , IN 
X( II )=YX( II ,  1  ) 

Y ( II )=YX( II , ICOL) 

8010  CONTINUE 

REWIND  INN 
C 

C  TERMINATE  DATA  READING  &  PLOT  CURVE 

C 

CALL  CURVE( X , Y , IN , 0 ) 

C 

C  RESET  LINE  TYPE 

C 

IF( ICT . EQ . 2  )  CALL  RESET(  ’ DOT ’  ) 

IF( ICT . EQ . 5  )  CALL  RESET(  ' DASH ’  ) 

IF( ICT . EQ . 4 )  CALL  RESET (’ CHNDOT '  ) 

IF( ICT . EQ . 5  )  CALL  RESET (’ CHNDSH 1  ) 

CALL  EXCMS( 'CLRSCRN'  ) 

WRITE( 4 , 9473  ) 

9473  FORMAT ( '  DO  YOU  WISH  TO  TERMINATE  THIS  SUBPLOT?’) 
READ( 4,9400)  ANS 
IF(ANS.NE.YES)  GO  TO  2000 
C 

C  END  OF  PLOT 

C  WRITE  ADDITIONAL  LABELS 

C 

CALL  ENDGR( 0 ) 

CALL  AREA2D( AX, AY) 

CALL  EXCMS(  ’CLRSCRN’  ) 

WRITE( 4 , 9480 ) 

9480  FORMAT ( ’  DO  YOU  WANT  ADDITIONAL  LABELS  ?',/) 

READ( 4 , 9400 )  ANS 

IF( ANS . EQ . ’ N ’ )  GO  TO  9900 

WRITE( 4 , 9483 ) 

9483  FORMAT ( '  YOU  HAVE  A  CHOICE  OF  FUNCTION 


+  -  1 - GAUSSIAN  FUNCTION’,/, 

+  ’  2 -  LEVEL  FUNCTION’,/, 

+  ’  6 -  TRAPEZOIDAL  FUNCTION’,/, 


+  ’  INPUT  THE  NUMBER  OF  FUNCTION  TYPE  (II)  :’,/ 

READ  (4,9401 )  NF 
IF  (NF.EQ. 1 )  THEN 

INF=  'GAUSSIAN  FUNCTIONS’ 

WRITE  (4,9484) 


FORMAT( 'SPECIFY  BUNCH  LENGTH  (F4.2 


-  1  0.24  CM  ’  ,/. 

■  2 - 0.88  CM 

■  3 - 1.20  CM  ’  ./. 

•  4 - 1.76  CM  ■  ,  /  , 

•  5 - 2.40  CM  *  ,  /  , 

■  6 - 2.63  CM  '  ,  /  , 

•  7 - 3.50  CM  •  ./. 

•  8 - 4.38  CM  ’  , 


'INPUT  THE  NUMBER ( II ) : _ ',/) 

READ( 4 , 940 1 )  NBL 
IF  (NBL.EQ. 1 )  THEN 

IBL=  ’BL  =  0.24  +CM ! $ ' 

ELSEIF( NBL.EQ. 2)  THEN 

IBL=  ’BL  =  0.88  +CM ! $ ' 

ELSEIF( NBL.EQ. 3)  THEN 

IBL=  'BL  =  1.20  +CM ! $ ' 

ELSEIF( NBL.EQ. 4)  THEN 

IBL=  ’BL  =  1.76  +CM ! $ ' 

ELSEIF( NBL.EQ. 5)  THEN 

IBL=  'BL  =  2.40  +CM ! $ ' 

ELSEIF( NBL.EQ. 6)  THEN 

IBL=  ’BL  =  2.63  +CM! $ ’ 

ELSEIF( NBL.EQ. 7)  THEN 

IBL=  ’BL  =  3.50  +CM! $ ' 

ELSE 

IBL=  'BL  =  4.38  +CM! $ ’ 

ENDIF 

ELSEIF( NF. EQ . 2 )  THEN 
INF=  'LEVEL  FUNCTIONS’ 

WRITE( 4 , 9486  ) 

FORMAT( 'SPECIFY  LEVEL  FUNCTION  LENGTH  (F4.2 ):’,/ 


>  1  - 1.00  CM  ’  ,  /  , 

■  2 - 3-00  CM  ’  ,  /  , 

•  3 - 4.50  CM  ’  ,/, 

■  4 - 5.00  CM  ’  ,  /  , 

'  5 - 6.00  CM  ’  ,  /  , 

’  6 - 7.00  CM  ’  ,/, 

>  7 - 8.00  CM  ’  ,/, 

’  8 - 9-00  CM  ’  ,/, 

'INPUT  THE  NUMBER(II): _ ',/) 

READ( 4 , 9401  )  LF 
IF  (LF.EQ. 1 )  THEN 

LFL=  ’ LFL  =  1.00  +CM ! $ ’ 

ELSEIF( LF . EQ . 2 )  THEN 

LFL=  ’LFL  =3.00  +CM ! $ ’ 

ELSEIF( LF . EQ . 3 )  THEN 

LFL=  ’LFL  =  4.50  +CM!$’ 

ELSEIF( LF.EQ. 4)  THEN 

LFL=  ’LFL  =  5-00  +CM 1  $ ’ 

ELSEIF( LF . EQ . 5 )  THEN 

LFL=  'LFL  =  6.00  +CM'$’ 
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ELSEIF( LF . EQ . 6 )  THEN 

LFL=  ’ LFL  =7.00  +CM1S’ 
ELSEIF( LF . EQ . 7 )  THEN 

LFL=  'LFL  =  8.00  +CM!$’ 
ELSE 

LFL=  'LFL  =  9.00  +CM1S' 
ENDIF 

ELSEIF( NF. EQ . 6 )  THEN 

INF  =  'TRAPEZOIDAL  FUNCTIONS’ 
WRITER, 9501  ) 

FORMAT ( 'SPECIFY 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


9502 


READ (4 


TOP 
1  -■ 


2 

3 

4 

5 

6 
7 


LEVEL  LENGTH  :  ’ 
--  0.07  CM  * ,  /  . 

--  2.94  CM  ' ,/, 

--  3.185  CM  ’ ,/ 
--  4.0  CM  ' ,/, 
--  5.0  CM  ' ,/, 
--  5.04  CM  ’ , 

--  8.96  CM  ’ 

NUMBER  ( II  ) : 


/, 


INPUT  THE 
,9401 )  TLL 
IF  (TLL.EQ. 1 )  THEN 

ITLL=  ' TFTL  =  0.07  +CM!$’ 
ELSEIF  (TLL.EQ. 2) 

ITLL=  'TFTL  =  2 
ELSEIF  (TLL.EQ. 3) 

ITLL=  'TFTL  =  3 
ELSEIF( TLL.EQ. 4) 


./) 


THEN 

.94  +CM ! $ ' 
THEN 

.185  +CM ! $ ' 
THEN 


ITLL=  'TFTL  =4.0  +CM! $ ’ 
ELSEIF( TLL.EQ. 5)  THEN 
ITLL=  'TFTL  =5.0  +CM!$' 
ELSEIF( TLL.EQ. 6)  THEN 


ITLL  = 

ELSE 
ITLL  = 

ENDIF 

WRITE( 4 , 9502 ) 
FORMAT ( 'SPECIFY 


TFTL  =  5.04  +CM ! $ ' 


TFTL  =  8.96  +CM! $ ' 


READ( 4 


BASE  LENGTH  : ' , / , 

1  - 5.215  CM  '  ,/ 

2  - 7.0  CM 

3  - 8.0  CM 

4  - 8.96  CM  '  ,/, 

INPUT  THE  NUMBER  (II): _ 

9401  )  BL 
IF  (BL.EQ. 1 )  THEN 

IBL=  ’ TFBL  =  5.215  +CM!$' 

ELSEIF( BL . EQ . 2 )  THEN 

I BL=  'TFBL  =  7.0  +CM 1  $ ' 

ELSEIF( BL . EQ . 3 )  THEN 

IBL=  'TFBL  =  8.0  +CM!$' 

ELSE 

IBL=  'TFTL  =  8.96  +CM!$’ 

ENDIF 


./) 
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'.V 
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ENDIF 

WRITE( 4 , 9485 ) 

FORMAT( ’  SPECIFY  RADIATION  MODE  (II):’.//, 

+  *  1 - CERENKOV ’ , /  , 

+  ’  2  SUB-CERENKOV’,/, 

+  ’  INPUT  THE  NUMBER: ’,/) 

READ( 4 , 940 1 )  NMOD 
IF(NMOD.EQ. 1 )  THEN 

IMOD  =  ’C+ERENKOV! $ ' 

ELSE 

IMOD  =  ’ +SUB-! C+ERENKOV! $• 

ENDIF 

WRITE(  4 , 94-87  ) 

FORMAT ( ’  SPECIFY  ENERGY  OR  NBETA  (II):’,// 
X  ’  FOR  E  ’  ,/, 

X  *  1  —  E  =  2.555  MEV* ./. 

X  ’  FOR  NBETA ' . / , 


FORMAT ( 


NBETA 

NBETA 

NBETA 

NBETA 

NBETA 


X  ’  INPUT  THE  NUMBER :_ 

READ(  4 , 94-01  )  NB 
IF(NB.EQ. 1 )  THEN 

INB  =  ’+N!*B!  =  1 .00$’ 
INB  =  'E  =  2.555  M+E ! V$ 
ELSEIF( NB . EQ . 2 )  THEN 

INB  =  ’+N!*B!  =  0.99$’ 
ELSEIF( NB .EQ . 3 )  THEN 

INB  =  ' +N ! %B !  =  0.98$’ 
ELSEIF( NB . EQ . 4 )  THEN 

INB  =  1 +N ! %B !  =  0.97$’ 
ELSEIF( NB . EQ . 5 )  THEN 

INB  =  ’+N!*B!  =  0.96$’ 

ELSE 

INB  =  ’+N!*B!  =  0.95$' 
ENDIF 

CALL  EXCMS( ’ CLRSCRN ' ) 
WRITE(4.9488) 


-.99  NM* ./, 
0.98  NM’ ,/ , 
0.97  NM’ . / , 
0.96  NM' ,/ , 
0.95  NM’ ,/ , 
:  ’  ./) 


FORMAT ( 


SPECIFY  ETA  ( II ) : 

1  --  ETA  = 

2  --  ETA  = 


ETA 

=  12 

ETA 

=  24 

’  ./. 

ETA 

=  57 

ETA 

=  74 

*  ./. 

<  ’  INPUT  THE  NUMBER: 

READ(  4 , 94-01  )  NETA 
IF( NETA . EQ . 1 )  THEN 
IETA  =  ' %C !  =  12  $’ 

ELSEIF( NETA . EQ . 2 )  THEN 
IETA  =  ’ %C !  =  24  $' 

ELSEIF( NETA . EQ . 3  )  THEN 


•  f  m  *  0  * 


IETA  =  ’*C!  =  37  $' 

ELSE 

IETA  =  ’itC!  =  74  S’ 

END  IF 

WRITE( 4 , 9490 ) 

FORMAT( ’  SPECIFY  BEAM  PATH  LENGTH  (II  ):*,//, 
C  ’  1  --  L  =  21  CM  (ETA=1 2 ,N=6 ) 1 

C  1  2  --  L  =  42  CM  (  ETA=  1  2  ,  N=3  ) 1 

C  ’  3  --  L  =  64.8  CM  (ETA=37,N=< 

(  ’  4  --  L  =  126.2  CM  ( ETA= 1  2 ,  N  = 

(  '  5  --  L  =  129.5  CM  (ETA=37,N= 

C  ’  INPUT  THE  NUMBER: _  ’,/) 

READ( 4 , 9401 )  NL 
IF(  NL  .  EQ  .  1  )  THEN 

INL  =  'L  =  21  +CM ! $ 1 
ELSEIF( NL . EQ . 2 )  THEN 
INL  =  ’L  =  42  +CM ! $ 1 
ELSEIF( NL . EQ . 3 )  THEN 
INL  =  ’L  =  64.8  +CM ! $ 1 
ELSEIF( NL . EQ . 4 )  THEN 


*  ./. 

6)’ ./. 
=1 )* ./, 
=  3)’  ,/. 


INL  = 

'L  =  126.2  +CM ! S ' 

ELSE 

INL  = 

’L  =  129.5  +CM ! $ ’ 

ENDIF 

WRITE( 4 , 

9492) 

format; ’ 

SPECIFY  WAVE  LENGTH 

(11  ): 

'.II, 

X  , 

1  --  RAMDA  = 

10.52 

CM  (N=1  )  1  ,  / 

X 

2  --  RAMDA  = 

3.50 

CM  ( N  =  3  ) ' , /  , 

X 

3  --  RAMDA  = 

1  .75 

CM  ( N  =  6  ) 1  ,  /  , 

(  ’  INPUT  THE  NUMBER 

READ( 4 , 9401 )  NRAM 
IF( NRAM . EQ . 1 )  THEN 
IRAM  =  ’*L!  =  10.52  +CM!  $ 1 
ELSEIF( NRAM. EQ . 2 )  THEN 
IRAM  =  ’ %L !  =  3.50  +CM ! $  1 


ELSE 
IRAM  = 
END  IF 


=1.75  +CM! $ ' 


XMES  =  2.880 
YMES  =  6.130 
I MESS  =  I MOD 

CALL  MESSAG( IMESS, 1 00, XMES, YMES) 


YMES  =  YMES  -  0.250 
IMESS  =  INB 

CALL  MESSAG( IMESS, 100, XMES, YMES) 


C 


YMES  =  YMES  -  0.250 
IMESS  =  IETA 

CALL  MESSAG( IMESS, 100, XMES, YMES) 


YMES  =  YMES  -  0.250 
IMESS  =  INL 

CALL  MESSAG( IMESS, 1 00 , XMES , YMES) 


YMES  =  YMES  -  0.250 
IMESS  =  IRAM 

CALL  MESSAG( IMESS, 1 00, XMES, YMES) 

C 

YMES  =  YMES  -  0.250 
IF(NF.EQ. 1 )  THEN 
IMESS= IBL 

ELSEIF( NF . EQ . 2 )  THEN 
IMESS=LFL 

ELSEIF( NF . EQ . 6 )  THEN 
IMESS=  ITLL 

CALL  MESSAG( IMESS, 1 00, XMES, YMES) 

YMES  =  YMES-0.250 
IMESS  =  IBL 
ENDIF 

CALL  MESSAG( IMESS, 100, XMES, YMES) 

C 

IF( HH . EQ . 6 )  THEN 
XMES  =1.20 
ELSE 

XMES  =1.50 
ENDIF 

YMES  =  -2.25 
IMESS  =  INF 

CALL  MESSAG( IMESS, 100, XMES, YMES) 

C 

9900  CONTINUE 

CALL  EXCMS( ' CLRSCRN ' ) 

IF( IOUT . EQ . 2 )  THEN 
WRITE(  4 , 94-93  ) 

9493  FORMAT ( ’  < CLEAR > ,  THEN  READY  FOR  HARD  COPY’) 
ENDIF 

CALL  ENDPL( 0 ) 

WRITE( 4 , 9495 ) 

9495  FORMAT ( ’  DO  YOU  WANT  TO  MAKE  ANOTHER  FIGURE?’,/) 
READ( 4 , 9400 )  ANS 
IF( ANS.NE.YES)  THEN 
GO  TO  9999 
ELSE 
ENDIF 

CALL  RESET (  ’ALL’  ) 

GO  TO  100 

9999  CONTINUE 

CALL  DONEPL 

STOP 

END 
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